Controlling the quantum state of a single photon emitted from a single polariton 
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We investigate in detail the optimal conditions for a high fidelity transfer from a single polariton 
state to a single photon state and subsequent homodyne detection of the single photon. We assume 
that, using various possible techniques, the single polariton has initially been stored as a spin-wave 
grating in a cloud of cold atoms inside a low-finesse cavity. This state is then transferred to a single 
photon optical pulse using an auxiliary beam. We optimize the retrieval efficiency and determine 
the mode of the local oscillator that maximizes the homodyne efficiency of such a photon. We find 
that both efficiencies can have values close to one in e 



I. INTRODUCTION 

The generation and characterization of well-controlled 
single photon states have accomplished considerable 
progress during recent years. Beyond the usual non- 
classical effects obtained using single photon counting 
and intensity correlation measurements [l[, interference 
effects offer very efficient ways to get more information 
about the quality of single-photon wave packets. For in- 
stance, observing coalescence, that is the famous "Hong- 
Ou-Mandel" effect @ , between two single photons emit- 
ted by two different atomic sources is a very good way to 
ensure that these two photons are indistinguishable Q. 

Another avenue in order to fully characterize a sin- 
gle photon state is using homodyne detection, and ap- 
plying quantum tomography techniques in order to re- 
construct the Wigner function W(x,p) of the photon in 
phase space, where x and p are the quadrature operators 
of the quantized electric field [|| . Such a method provides 
a full characterization of the quantum state of the sin- 
gle photon, which is very intuitive because the quality of 
the single photon state, and especially its purity, directly 
translate into the negativity of the Wigner function. 

In the optical domain quantum homodyne tomogra- 
phy has been implemented for various number states, in- 
cluding one- 3 and two- [!, Q photon Fock states. In 
these experiments, the usual method to produce a single- 
photon state is to use parametric fluorescence emitting 
non-degenerate pairs of photons, and to detect one of 
them using a photon counter, such as an avalanche pho- 
todiode. Then the other beam is projected in the desired 
one-photon Fock state, which is analyzed by homodyne 
detection and then quantum tomography. It has to be 
noted that the quantum efficiency of the homodyne de- 
tection channel has to be very high: any loss at that stage 
degrades the purity of the single photon state by mixing 
it with vacuum, and quickly destroy the negativity of the 
Wigner function. This method works quite well, but it 
is intrinsically non-deterministic : the probability of the 
first (photon counting) event must be low, and then the 
photon cannot be emitted "on demand" , when needed 
for applications. 

Here we would like to investigate another scheme, 



FIG. 1: (Color online) An ensemble of N three-level atoms 
placed in a cavity is initially prepared in a single polariton 
state, with 1 atom in the metastable state |s) and (N — 1) 
atoms in the ground state \g). The atoms are coupled to the 
quantum field with collectively enhanced coupling w = gvN, 
where g is the single photon Rabi frequency. By applying a 
classical control field (red) with Rabi frequency Q(t) a single- 
photon state is generated with well defined properties. The 
cavity detuning S is used to optimize the photon retrieval 
efficiency. 



which can, at least in principle, be made more deter- 
ministic. The idea is to emit the single photon from a 
so-called "polariton" state, where a single long-lived spin 
excitation is distributed over many atoms. It can be con- 
verted into a free-propagating single photon using a con- 
trol laser beam (see Fig.[T]). As we will show below, given 
the presence of the polariton the source is deterministic: 
the phase-matching condition between the readout laser 
beam, the atomic spin wave and the single photon leads 
to a collective enhancement effect which ensures that the 
photon is emitted in a well-defined spatial and temporal 
mode. For the whole scheme to be deterministic, the po- 
lariton should also be prepared in a deterministic way, 
which can be done for instance by using the so-called 
"Rydberg blockade" mechanism as explained in recent 
theoretical papers @, Q. However, even when the prepa- 
ration of the polariton is probabilistic, one may exploit 
the fact that the single photon can be emitted on de- 
mand, just when it is actually needed, within the rela- 
tively long coherence time of the spin wave. In this article 
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we will not consider in more details how the polariton is 
prepared [9l-fl3j. but rather look at the following ques- 
tion: assuming that a single polariton is imprinted on 
an atomic cloud within an optical cavity, how efficiently 
can we turn it back into a single photon within a single 
mode? What will be the amplitude and frequency of the 
emitted single photon wave-packet, and is it possible to 
control it, and to mode-match it efficiently onto a local 
oscillator? In order to answer these questions, we will use 
as a starting point the theoretical analysis presented in 
which has been shown to be in good agreement with 
experiments [l4j . We will extend this analysis in order to 
characterize very carefully the homodyne efficiency which 
can be obtained for such a single-photon pulse, since this 
efficiency is crucial for quantum homodyne tomography. 
Finally we will discuss various experimental and theoret- 
ical considerations. 



II. THEORY 

We assume that the atoms are confined in an optical 
cavity and a single polariton state has been prepared so 
we only analyze the retrieval of a single photon. We are 
primarily interested in situations in which the collectively 
enhanced coupling between the atoms and a quantized 
cavity radiation mode is always one of the dominant en- 
ergy scales, but not necessarily the only one. This is an 
interesting regime for experiments and it has not been 
studied in detail before. The underlying atom-cavity 
model has been already described @ and thus we only 
briefly introduce it. 



A. Cavity model 

Various techniques and strategies for photon storage 
and retrieval have been analyzed using a relatively simple 
model [(|. It is assumed that a cloud of N three- level 
atoms (see Fig. [IJ with two metastable states \g) and 
\s), and an excited state |e) are coupled to a classical 
and quantized single-mode fields. The classical (read) 
field couples the states |s) and |e) while the quantized 
field couples the states \g) and |e). The coupling to the 
classical field is characterized by a slowly varying Rabi 
frequency f2(t) and the collectively enhanced coupling w 
between the atoms and the quantized cavity radiation 
mode is w = gV~N, where g is the single photon Rabi 
frequency. 

The assumption that, initially, a single excitation is 
shared by a very large number of atoms practically 
means that almost all atoms are in the ground state 
at all times so that the following approximate relations 
are justified: » N, and Z? = i<e « ££i<4 « 

^^=i a ls ~ 0, where <j l uv = \u) ti (v\. Under these 
simplifications and in the rotating-wave approximation, 
we define the atomic polarization annihilation operator 



P(t) = '^2^Li fge/V^V", and spin-wave annihilation opera- 
tor S(t) = J2iLi a% gs l V^V- In the general case, these oper- 
ators also include phase factors exp(— i k.n), where k is 
the laser wave-vector, and ri the position of atom num- 
ber i. There phase factors can be absorbed by redefining 
the state \g)i- In Appendix A we write the Langevin 
equations for the coupled propagation of the field annihi- 
lation operator £(t) and atomic operators P(t) and S(t) 
in the Heisenberg picture, and we show that it is actually 



enough to consider the c-numbcr equations : 

t = -(n + iS) £ + iwP, (1) 

P = -('y + iA)P + iw£ + iflS, (2) 

S = iWP, (3) 



where k and 7 are respectively the cavity and polariza- 
tion decay rates, while 5 is the detuning between the 
laser-driven Raman light and the cavity, whereas A the 
detuning between the driving laser and the atomic line 
(see Fig. [J). 

It should be pointed out that the c-numbers £(t), P(t) 
and S(t) are not averages of the corresponding operators, 
but can be seen as quantum amplitudes, which allow one 
to evaluate normally ordered products of operators (see 
details in Appendix A). An intuitive justification for this 
approach is that, since there is only one excitation in 
the system, there will be only one emitted photon in the 
cavity mode; we are therefore calculating the amplitude 
of this single photon wave packet. The initial conditions 
are £(io) = P(to) = and S(to) = 1, where to indicates 
the beginning of the retrieval sequence. 

The retrieval efficiency rj is defined as the ratio of re- 
trieved photons and stored excitations. Since the number 
of stored excitations is assumed to be one, 77 is equal to 
the number of retrieved photons 

/>oo 

7/= / dt\£ out {t)\ 2 , (4) 

Jto 

where 

£ out (t) = V2^£(t). (5) 

As said before, we do not need to use quantum operators 
in Eq. Q since we never have more than one photon in 
the quantized radiation mode. 

The propagation equations for (£r,Pi,Sr) and 
(£1 , Pr, Si) decouple for A = 5 = and real 0, where 
(1r and ai are the real and imaginary part of a com- 
plex variable a. Consequently, £ and S are real and P is 
imaginary in that case. For real H, there is the following 
mapping between the solutions for ±A 

£(-A,-5,t) = £*(A,5,t), (6) 
P(-A,-S,t) = -P*(A,M), (7) 
S(-A,-5,t) = S* (A, S,t). (8) 

For non-zero A, the interaction between the atoms and 
the quantized radiation mode shifts the cavity resonance. 
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Therefore, to compensate this shift and maximize the ef- 
ficiency Q, we introduce some cavity detuning 5 7^ 0. 
Relations ([6])- ((8]) impose the mapping between the opti- 
mized cavity detunings S op t as follows 



S opt (-A) = -<5 opt (A). 



(9) 



As explained in Section [HI! other mappings relevant for 
the optimization of homodyne detection (HD) also follow 
from Eqs. All these mappings also have practical 

usage since they reduce numerical work by half. 



B. Retrieval efficiency 

We want to explore the parameter space by means of 
numerical calculations to find regions of high photon re- 
trieval efficiency. In many experiments the emission of 
the single photon is detected by a counter, such as an 
avalanche photodiode, and the relevant parameter is i] 
as defined by Eqs. [4] and [5l Here, we arc also inter- 
ested in using homodyne detection, which allows us to 
analyze fully the single photon wavepacket, and e.g. to 
reconstruct its Wigner function by using quantum opti- 
cal tomography. Therefore we will also look for optimal 
conditions for homodyne detection. 

Let us first review the expression for the maximal pho- 
ton retrieval efficiency in a general case since our search 
for optimal retrieval conditions do not necessarily involve 
any special cases such as the bad-cavity limit (n ^> w). 
Actually, based on experimental considerations, the typ- 
ical range of parameters we are interested in corresponds 
to if >> k ~ 7. We will show now that the efficiency 
maximum is then the same one derived in @, though we 
are no more in the bad cavity limit. 

A common way to simplify the problem is to try to 
find the first integrals of the equations of motion. Since 
our system is dissipative (k, 7 ^ 0), we may at least look 
for some relations analogous to the first integrals. The 
following relation gives just that 



l^| 2 + |S| 2 ) =-2 7 |P| 2 -2 K |£f 



(10) 



We are interested in the conditions for high retrieval 
efficiencies r\. Using Eqs. ([5]) and (fTU|) . as well as the 
conditions S(t ) = 1, £(t ) = £(00) = P(t ) = P(co) = 
5(00) = 0, one gets the following relation 



1 = 77 + 27 / \P(t)\ 2 dt. 



(11) 



The boundary condition S (00) = is only fulfilled for 
maximal efficiencies since the spin- wave decay rate is ne- 
glected. It is convenient to rewrite the integral in the 
last equation in the frequency space and then use the re- 
lation between the Fourier components P(lu) and £(u>) 
obtained from Eq. (TTJ) 



w 



(12) 



Therefore 

\P{t)\ 2 dt 



\P(co)\ 2 duj 



to 



(co-S) 2 



|£(w)| 2 dw,(13) 



which implies 



2 7 / \P(t)\ 2 dt^ 



2k / \£{t)\ 2 dt 



C 



(14) 



where C = ui 2 /k7 is the coopcrativity factor. Combining 
this condition and Eq. (fTTj). we get the efficiency limit 

«*TTc- (15) 

As shown in Q for 5 = 0, it is possible to achieve the 
highest efficiency C/(l + C) in some particular regimes. 
Since the relation (|T5|) is very general, introducing an op- 
timized cavity detuning S opt ^ may allow us to optimize 
the efficiency, especially for low intensities of the driving 
field, but it will never exceed this limit. 

A necessary, but not sufficient, condition to have a 
maximal efficiency is that the term (uj — S) 2 /k 2 in Eq. 
(fT3"|) should be negligible. For 6 = 0, this happens if the 
range of frequencies uj corresponding to non-vanishing 
P(uj) and £(uj) is much smaller than k. In this case, 
achieving high efficiencies would require long read pulses 
or high Q. However, according to Eq. (fl"3|) . the efficiency 
can be improved significantly by introducing a cavity de- 
tuning 5 that cancels uj corresponding to the maximum 
of A consequence is that S op t is comparable with 

the spectral width of |P(w)| and so that we expect 

that both uj 2 /k 2 and 8 2 pt /n 2 are negligible for very long 
pulses. Because the term uj 2 /k 2 in Eq. (| 13[) originates 
from £(t) in Eq. JT]), we may think that £{t) can be 
neglected for very high efficiencies. Dropping this term 
yields to the 'bad cavity' relation between £{t) and P(t) 



— k£ + iwP = 0, 
and therefore (k, w > 0) 

k\£\ = w\P\. 



(16) 



(17) 



In general, as shown in detail in Appendix B, for suffi- 
ciently long pulses and high efficiencies a better approx- 
imation can be obtained by noting that Eq. (fTJ is the 
first-order expansion of 

- n£(t + St) + iwP(t) = 0, (18) 

with respect to the short constant time delay St = 1/k. 
Equation (fTS)) substitutes to Eq. ((T0|) in the time domain, 
while in the frequency domain it becomes 

- e i,p "K£(u}) + iwP{uj) = (19) 
with W = —uiSt, and the relationship (|17|) remains valid. 
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C. Solution for large w 

In this section, we assume that w is at least large but 
f2 and A can be large as well. For w 3> k, 6, relation 
(|T2"j) shows that the ratio |P|/|£| is small as long as the 
pulscwidths of ft, £, and P are of the same order of mag- 
nitude 



|P| 



|P| 



I -PI 

- — -«1 

\£\ « 



(20) 



Because S ~ 1 at short times (due to the initial condi- 
tion S(to) = 1), Eqs. ([20]) and © lead to the following 
estimate of £ for small A, 
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Let us introduce a new variable Q = £ + pP and take 
p for which the following equation is held 



Q = XQ + ipnS. 



(22) 



There are two combinations of p, X for which the last 
equation is fulfilled 



A.. 



7 + iA + fc± i v /4w 2 -(7+iA - k)' 



j+iA-k ± iy/4w 2 — (j+iA - k)' 



(23) 
(24) 



where k = n + iS. We choose p, X which facilitate the 
adiabatic elimination of Q in Eq. (f2"2"j) 



£ 



■ p p = -' l 2ns = ~ns, 



where 



X 

-iw— 
P 



(25) 



(26) 



is introduced for convenience. 

Estimates (|20p and (j2"Tj) show that |P| is a much 
smaller quantity compared to both \£\ and \S\. There- 
fore, dropping the small term pP in Eq. (|2"5)) has little 
effect to the relationship between much larger £ and S 



e—ffis. 



(27) 



The last equation can be obtained in a different way and 
somewhat different approximations. For some parame- 
ters, we can use both X± and p± to eliminate P in Eq. 
(f2"5l ) and get relation (|27p again with a new definition 

for /3 



-iw 



X+ - X- p+p- 



(28) 



For our typical range of parameters, there is no practical 
advantage of using either Eq. (|2"6"|) or Eq. (|2"5|) but it 
turns out that the definition (|28[) is more suitable for very 



FIG. 2: (Color online) Time delay between P{t) and £(t). (a) 
For high retrieval efficiencies and small Cl/Sl (i.e long pulses), 
the cavity field £ is delayed with respect to P by a constant 
delay 8t — 1/k: £ (t) ~ iwP(t — St)/n. The numerical solution 
for P (solid red) is compared with Pbd(i) = —in£(t)/w (solid 
green) and P B hat(t) = Pbd(t + St) (dashed blue), (b) For 
shorter pulses (r = 150 ns) this simple delay approximation 
breaks and P(t) and £ (t) acquire different shapes. 



long pulscwidths, as shown on Fig. [5J Since S « 1 in the 
beginning of a pulse sequence, we see that, for initial 
times, the shape of \£\ and are always the same and 
\£\ <C |fi|. Taking the derivative of Eq. (|27|) with respect 
to t and using Eqs. {T]) and (j3]) to eliminate S and P, we 
get 



/3 2 



£ 



£. 



(29) 



For real f2, the solution of this equation consistent with 
the initial condition £(irj) = is relatively simple 



/3 ^W+n- 



: exp 



r> 2 



/* 2 



n 2 



rdf' 



(30) 



where the multiplicative constant {—w/0) has been de- 
termined by taking limits of Eqs. (f3"0")l and (|27l) at t — >• to 
and utilizing the initial condition 5 (to) = 1- Similarly, 
the solution for S is 



S(f) 







: exp 



-k 



n 2 



n 2 



VW+W L J to 

Finally, the solution for P in our approximation is 

n 2 



i8(kn + n) 



(/? 2 + n 2 ) 



-k 



t 



2 + n 2 



-d£ 



(31) 



(32) 



Let us notice that if |fJ| 2 < |^| 2 , then w 
—wQ,(t)/f3 2 = —ipCl(t)/X, which means that ft and £ 
have the same pulse shape. However, according to Eq. 
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FIG. 3: (Color online) Comparison between numerical (solid 
red) and analytical (dashed blue) solutions for the emitted 
signal field |£(£)|. The Rabi frequency Q(t) of the driving 
field is also shown (dashed green), scaled down to match the 
pulse area of the corresponding cavity field |£(t)|. Let us 
note that the (very small) emitted pulse lies entirely within 
the envelope of the (very large) driving pulse, but its peak 
value does occur before the driving field peak value. The 
parameters are Q./2-k = 80 MHz, C = 200, 8 = MHz, and 
pulsewidth (FWHM) r = 200 ns for several detunings. The 
agreement is practically exact for A = 0, but it gradually 
worsens as A increases. 



P?]) and (j3"Tj) . this also means that S ps 1 and hence the 
retrieval efficiency is very low. 

Solutions ([3T) |) - (j3"!?)) are excellent approximations for 
large w and A sa 0, they work rather well for moderate to 
large A, but for large detuning they become inaccurate, 
as illustrated in Fig. [31 The reason why the approxima- 
tion worsens depends on which definition of (3 is used. 
For A = S = 0, both |A±| are the same but they can be 
very different for larger A. The same statement holds for 
|p±|. Since the larger \p\ corresponds to the larger |A|, 
dropping the term pP in Eq. (fz?5j) becomes problematic 
for larger A even though the adiabatic elimination itself 
works better. On the other hand, the problem with using 
Eq. (|28p for larger A is that the adiabatic elimination 
corresponding to the smaller |A| becomes unreliable. 



tocurrent i(t) is proportional to the following operator 

i(t) oc (a* a{t) + a a\t)) /\a\ = e ie a{t)+ e -* s a^t) (33) 

which is a quadrature operator with respect to the signal 
field, determined by the choice of 6. 

In the present work, we will consider a situation with 
a continuous single- mode LO, but where the signal beam 
corresponds to a single photon, which is an essentially 
transient field, triggered by the "read-out" pulse sent 
onto the atoms in the cavity. This situation is somehow 
similar to the homodyne detection of transient homodync 
signals followin g a p hotodetection event, as measured for 
instance in [l6|, 17 1. The idea is then to "extract" from 
the stationnary noise coming out from the homodyne de- 
tection a non-stationnary signal, through an appropriate 
time-domain filtering, which has to be matched as well 
as possible to the expected "single-photon wavepackct" . 

Two effects can essentially degrade this matching : the 
first one is the efficiency i] introduced in the previous 
section, i.e. the probability that the photon is effectively 
emitted in the output mode a(t), and the second one 
is the quality of the interference between the photon and 
the LO, represented by an overlap factor \ defined below. 

In order to examine the effect of the quality of the in- 
terference, we will neglect all technical mode-matching 
issues, and consider only the effect of the interference 
of the LO with the single-photon amplitude £(t) intro- 
duced in the previous section. Our goal is to optimally 
match the mode £det(t) measured by the homodyne de- 
tector to the mode of the single photon defined by £{t). 
The amplitude f(t) = \£det(t)\ can be controlled either 
by modulating the intensity of the local oscillator or by 
rapidly sampling the signal and applying a numerical fil- 
ter f{t) to the data, and in the following we will assume 
that the amplitudes of the two modes are identical: 



/(*) 



|£out(*)| 

Vv ' 



(34) 



III. OPTIMIZATION OF HOMODYNE 
DETECTION 

Homodyning consists in interfering two light fields: 
the local oscillator (LO), and the measured "quantum" 
beam (or signal beam), associated with mode annihila- 
tion operator a(t). The local oscillator is assumed to be 
monochromatic and to be in a coherent state \a), with 



I ale where 



is the LO photon flux and 9 is 



a controllable LO phase. More precisely, the signal from 
the homodyne detection (HD) is obtained by taking the 
difference of the two output photocurrcnts, resulting from 
the interference of the signal and LO. This difference pho- 



The phase 6(t) of £det{t) includes a constant phase shift 
ipo, which can be introduced on the local oscillator using a 
piezoelectric actuator, and correspond to various choices 
for the measured quadrature. In the present situation 
of a single photon state, we don't expect this phase to 
have any effect on the measured signal. In general a 
rotating phase ui^t can be added by changing the LO 
frequency by uiq. More complex functions 9(t) can be 
realized in principle, for example by chirping the local 
oscillator frequency, but for simplicity we will assume 
that 6(t) = tpo + ujot. 

The mode-matching efficiency can be optimized by 
maximizing the variance of the homodyne detection out- 
put, integrated over the output pulse. Up to a multiplica- 



G 



tive constant this variance is given by the expression 



V 



dtf(t)i(t) 



(35) 



= j dtdt'f(t)f(t') (i(t)i(0) 
= J dtdt'f(t)f(t')x 

((a(t)tf(t')) e id ^- i9 ^ + (ot(t)o(f )) e id ^- id W + 
{a{t)a{t')) c^+^t") + (tf{t)tf(t')) e -m-m')) 

In the present case, the initial conditions are such that 
the output cavity field will be a diagonal mixture of 
zero and one photon, and therefore the average values 
(a(t)a(t')) and (a^(t)d^(t')) are zero for all times (see 
discussion in Appendix A). Since the calculated quan- 
tities provide normally ordered averages, one uses the 
commutation relations to get finally : 



V = 1 
= I 

where 



2 J dtdt'f{t)f{t'){a){t)a{t'))e ie{t '^- ie{ -im) 
2I 



dtf{t)£ mt {t)e- iUot 



(37) 



As expected V = 1 (vacuum noise) if Jo = 0, and V = 3 
(single photon state) if Jo = 1. It should be noticed that 
the phase ipo does not appear any more (as it could be 
expected), and one has thus only to optimize V with re- 
spect to ujq. The optimal uiq, denoted by w pt, is obtained 
by maximizing the integral in and is found numerically. 
Up to a constant prcfactor, the optimized detector signal 
7 max is equal to the maximal value of Iq. Interpreting 
the integral in Eq. (|37[) as a scalar product of f{t)e lu ' ot 
and £out{t), and using the Cauchy-Schwarz inequality for 
real and normalized (||/|| =1) filters, we get 



< 



|£|| 2 = J dt\£ 



t (*)| 2 



(38) 



The upper bound / upP er = V of In, given by Eq. 
can be used to define a measure x °f t ne optimization 
success, by the expression 

max(/ ) = |J rft/(t)£o U tft)e-'"°^| 2 
V 



X = 



-'upper 



The mapping (JSJ implies the following mapping 
optimized frequency shift 



(39) 
of the 



^opt 



(-A) 



^opt 



(A), 



(40) 



Nearly maximally allowed value 77 of Iq is achieved by us- 
ing the filtering function f(t), and by looking for the best 
possible linear (in time) approximation of the phase of 
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FIG. 4: (Color online) Linearization of the phase de of the 
quantum field £(t). To get the highest homodyne efficiency, 
the phases of the quantum field and local oscillator have to 
be matched (see Section III). Assuming the phase #lo of the 
local oscillator to be linear in time, the optimal frequency 
td op t of the LO and the optimal cavity detuning 8 op t that 
maximize the homodyne efficiency XV are numerically found 
for the generic parameters r = 150 ns and C = 100 used in 
section [TV] In panels a) and b), 9$ (red solid) and c9lo (blue 
dashed) arc shown for A/2tt=120 MHz and f2/27r-=13 MHz. 
No cavity detuning optimization (S = 0) is performed in case 
a) while an optimized <5 op t/27r =15.5 MHz is used in case b). 
The homodyne efficiency \ starts to decrease for a given A 
and sufficiently large Q because 9g becomes more complicated, 
as shown in panel c) for Q,/2n = 80 MHz and A/2n = 300 
MHz, with <5 op t/27T =—2.8 MHz. Panel d) is the same as c), 
except for the cooperativity that becomes C-=200, improving 
the phase matching quality. 



£out(t) (see Fig. 2]). As shown below, the regions of high 
values of r\ and x have a very good overlap, which is im- 
portant information for future experiments. Let us em- 
phasize however that in practice, £ ut(i) in Eq. (|39p has 
to be corrected by including experimental losses, which 
will result in lower overall detection efficiencies. 



IV. RESULTS 

In this section we present numerical results on optimal 
conditions for efficient retrieval and subsequent detection 
of a single photon. The efficiencies 77 and x are evaluated 
using the numerical solution for £ (t), assuming a read- 
out pulse with a Gaussian shape (as a function of time) . 
Based on a typ ical experimental design using Rb atoms 
in a cavity [la [. we will take k/27t = 9 MHz and 7/27r = 3 
MHz as fixed values, and explore a parameter space 
consisting of (A, f2) points. Similarly, we will assume 
that the cooperativity factor is in the range C ~ 100 
lla l. using for instance the two hyperfme states \g) = 
--2,M F = 2) and \s) = |5s 1/2 , F= 1, M F = l) of 
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FIG. 5: (Color online) Retrieval efficiency r\ (panel a) and 
overall homodyne efficiency x 7 ! (panel b) for non-optimized 
cavity detuning (i.e. 5 = 0) for C = 100 and r = 150 ns. The 
frequency shift cjo is optimized to get the best XV f° r anv 
given A and Q. 



FIG. 6: (Color online) Retrieval efficiency r\ (panel a) and 
overall homodyne efficiency x 7 ! (panel b) for optimized cavity 
detuning. The efficiency XV is optimized with respect to fre- 
quency and cavity detuning (same parameters as in Fig. [5J. 



Rb. The readout pulse duration is chosen short enough 
(150 ns) so that the motional decoherence can be ne- 
glected for a cold atom cloud flil . flo| . 

A. Optimization of the cavity detuning 

As pointed out in Section III Bt if A 7^ the frequency 
of the cavity mode can be shifted due to the presence of 
atoms. To compensate for this shift, we introduce a cav- 
ity detuning S and numerically find S = S opt that maxi- 
mizes either the efficiency 77 or the product XH- The most 
dramatic improvements occur at lower fi, as shown by 
comparing Figs. [5]and[6j and the dependence 5 op t(A, il) 
is shown in Fig. [7] for C= 100. The optimum cavity de- 
tuning for large f2 is rather small, but it further flattens 
the high efficiency "plateau" visible on Figs. [5] and [6] 

Optimizing XV gives the upper homodyne detection 
limit, according to Eq. (j3"9")l . Ideally, the retrieval and 
detection efficiencies should be very close to each other 

max[x(S, u!o)rj(d)] = max [77(5)]. (41) 

Since the right hand-side of the last relation only depends 



on 5, we can first find <5 pt that optimizes rj(6) and then 
subsequently optimize x(<Wt 7 wo)ry(<5 pt) with respect to 
u>q. This optimization procedure can be done very ef- 
ficiently since only one parameter is varied at the time. 
This <5 op t corresponds to the maximum retrieval efficiency 
and thus it is an useful parameter in itself. The values 
5 opt and w opt obtained in this optimization step can then 
be used as an initial guess to find optimal S and ujq for XV- 
However this second optimization step does not lead to 
any practical increase in xVi as long as the linear phase 
approximation is valid (see Fig. |4|). In that case, the 
optimization of XV with respect to S and ujq can be effi- 
ciently replaced by optimizing ri(8) with respect to S and 
subsequently x((5 op t,wo) with respect to ujq- 

For example, in the case presented in Fig. (J4jb (resp. 
(|4])c), the x^-optimization gives (5 op t/27r = 15.5 MHz 
(resp. -2.8 MHz) while the ^-optimization gives <5 op t/27r 
=15.7 MHz (resp. -5.2 MHz). But even in the appar- 
ently less favourable case of Fig. Qc, the products XV 
are almost the same and nearly maximal, because the 
decrease in x is compensated by an increase in 77, so that 
XV = 0.930 (resp. x 7 ? — 0.937) as shown on Fig. SJ 
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FIG. 7: (Color online) Optimized cavity detuning for r=150 
ns and cooperativity C — 100, which corresponds to w/2n= 
46.5 MHz. The cavity detuning <5 opt that maximizes the effi- 
ciency xV is found for any given (A, fl). 



B. Calculation of the efficiencies 77 and xq 

Numerical results for 77 and x 7 ! are shown in Fig. [5] and 
|forr= 150 ns and C = 100. On Fig. one has S = 
(i.e. no (^-optimization is performed), whereas Fig. [5] 
corresponds to the optimized cavity detuning <5 opt (A, f2), 
shown on Fig. [Jj We observe that both quantities are 
symmetric under the transformation A — > — A, consis- 
tently with the mappings © and (J7J. 

Without cavity detuning optimization (see Fig. [3]), 
a plateau of nearly maximum retrieval efficiencies is 
reached for sufficiently large Rabi frequencies f2 of the 
read field. Increasing the pulse duration t helps estab- 
lishing the efficiency upper limit C/(C+1) in a larger 
region in the parameter space, because the spectral dis- 
tribution becomes more narrow for longer pulses so 
that ~ oj 2 /« 2 in Eq. (fl~3)l is negligible. Increasing the co- 
operativity factor increases the efficiency limit C/(C+1) 
but also widens the plateau of nearly maximal efficien- 
cies. Since we keep k and 7 fixed, larger C values imply 
larger couplings w. However, for A w 6 ~ 0, the growth 
rate of 77 as 51 increases is lower for larger w because, ac- 
cording to the solution ([50]) . £ is a function of the ratio 
51//3 and /3 ~ w. Consequently, we need to operate with 
more intense read fields for larger C in order to retrieve 
single photons with near certainty. For example, increas- 
ing C from 100 to 200 for A = 8 = requires ~ 40 % 
higher Rabi frequencies to reach the plateau region. 



C. Optimized values of the efficiencies 77 and \r\ 

In Fig. [6l we show ?7(A, 51) assuming that 5 is opti- 
mized. As expected, the cavity optimization yields to 
large improvements in efficiency for smaller 51. As Figs. 
[S^, and [5^ indicate, there is a region of relatively small A 




FIG. 8: (Color online) Frequency shifts = a> op t/27r corre- 
sponding to optimized cavity detuning (the curve is antisym- 
metric with respect to the line A = 0, v — 0). These values 
of v optimize the overall homodyne efficiency x 7 ! in a linear 
phase approximation (see Fig. [6}. The parameters are the 
same as in Fig. [7J 

and 51 in which the ^-optimization can significantly boost 
77 to reach almost 1, as shown for instance in Fig. [4]a,b. 

An overall effect of the (S-optimization to XV can be 
seen by comparing Figs. [SJd and [BJd. As a general rule, 
the curve on Figs, ^jp broadens and flattens, and high 
efficiency values become available for lower readout Rabi 
frequencies tt. In general, increasing any of C and r 
expands the region of maximal XV- Finally, Fig. [S] shows 
the behaviour of ^ opt = oj op t(A, D,)/2tt. The frequencies 
u op t are rather small and of the order of a few MHz. For 
A = and real J7, £{t) is real and negative and thus no 
phase correction is necessary This could be a preferable 
experimental choice, provided that a high enough Q is 
available. The function w op t changes sign across the line 
A = due to the symmetry f4"0")) . As previously pointed 
out, the function is opt exhibits distinguishable features 
for relatively small A and 51 due to the cavity detuning 
optimization. 

V. CONCLUSION 

As a conclusion, we have calculated the amplitude and 
phase of a single-photon wave-packet, obtained by de- 
tcrministically "reading out" a single polariton, which is 
assumed to be written beforehand onto a cloud of atoms 
within a low-finesse optical cavity. This "writing" can 
be done either non-deterministically by scattering a pho- 
ton from the atoms [h4T |. or deterministically by using 
Rydberg blockade techniques 0. The results indicate 
that both the emission efficiency and the homodyne ef- 
ficiency towards an adequate local oscillator can reach 
values very close to one. This opens the way towards the 
full quantum homodyne tomography of a single photon 
deterministically generated within a single mode. 
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Appendix A 

As shown in rcf. Q, the coupled propagation of the 
field annihilation operator £(t) and atomic operators 
P(t) and S(t) in the Heisenberg picture writes : 

£ = -(k + iS)£ + iwP + V2k £ in , (42) 

P = -(-/ + iA)P + iw£ + iQS+ y^Tp, (43) 

S = i£l*P, (44) 

where k and 7 are respectively the cavity and polariza- 
tion decay rates, while S is the detuning between the 
laser-driven Raman light and the cavity, and A the de- 
tuning between the driving laser and the atomic line (see 
Fig. [T]). Here we omit the spin- wave decay rate, but it 
can be reintroduced easily as shown in Q . The (bosonic 
vacuum) Langevin noise terms are £i n and Tp. 

In the article, we treated £ , P or S as complex numbers 
and ignored all noise terms. Here we want to explain 
what c-number quantity £ (t), P{t) and S(t) represent 
and in which situations they can(not) be used. First, 
let us note that £(t), P(t) and S(t) are not expectation 
values of the corresponding operators because 

(£{t)) = (P(t)) = (S(t)) = (45) 

is fulfilled at all times. The last expression is true because 
there is only one excitation in the system, which may 
be in either of the bosonic operators £ , P or S. Our 
main goal in this appendix is to show that £(t), P{t) 
and S(t) determine the expectation values of normally 
ordered operators as follows 

(Al(t)A j (t)) = A i (tyA j (t), W 

where Aij can be any of the operators £, F or S. We 
point out that the noise terms vanish in the propaga- 
tion equations for dt(A\(t)Aj(t)) because we assume that 
the input noise terms £ j„ and Tp are in the vacuum 
state. Then it is easy to check that all normally ordered 
noise correlations are zero @. This means that, as far 
as the averages (A\(t)Aj(t)) are concerned, the propaga- 
tion equations for dt{A\(t)Aj(t)) are the same ones as for 
d t (A*(t)A 3 (t)), which can be obtained from Eqs. (|42)l - 
(|4"4"|) with all noise terms set to zero. Using this conclu- 
sion, the relation (|46[) follows from the following simple, 
but nonetheless important mathematical statement: If 
the evolution of operators Ai and Bi , 1 < i < n is linear 

Ai = Y,^s(t)A s (t),, (47) 

s 

k = 53 (48) 



and if the expectation values of all products (BiAj) at 
some t = to can be factorized as 

(Bi(to)Aj(to)) = B m A j0 (49) 

then (Bi(i)Aj(t)) are factorisable at all times 

{Bi(t)M*)) =B i (t)A j (t), (50) 

where Bi(t) and Aj(t) satisfy Eqs. (j37J and (jUJ with 
the initial conditions Bi(to) = -Bio an d Ai(to) = AiQ. 

We can directly check that the products Bi(i)Aj(i) 
always satisfy the equation of motion for (Bi(t)Aj(t)) . 
However, to be solutions for (Bi(t)Aj(t)} , the products 
Bi(t)Aj (t) have to also satisfy the initial conditions. This 
is only possible if all (Bi(t )Aj(t )) are factorisable as re- 
quired by Eq. (|4^1) . We emphasize that this factorization 
of (Bi(to)Aj(to)) is not a trivial condition for operators 
and it is not generally fulfilled. 

In our case, the only initial nonzero expectation value 
is {§' (to) S (to)) so the initial averages are clearly fac- 
torisable. This means that the condition (|4"9"]l is fulfilled 
and so is the relation (|46)). Therefore, we have shown 
that the expectation values of the normally ordered op- 
erators are exactly the same as for the products of the 
c-number quantities £(t), P(t) and S(t). For simplicity, 
S(to) = 1 is used because the only effect of another phase 
convention S(t ) ->■ S(to)e lclH < is that all £(t), P(t), and 
S(t) acquire the same additional phase factor e 1 ^ which 
leaves the averages (|46p unchanged. 

In case where there is only one excitation in the sys- 
tem, the condition (|49|) is equivalent to the requirement 
that the initial state is a pure state. To show this, 
we consider here the evolution of the density matrix of 
the system (in the Schrodinger picture), rather than the 
evolution of the operators (in the Heisenberg picture) 
used up to now. We take as a basis the Fock states 
{|000) , 1 100) , |010) , |001)} for the bosonic operators £, 
P or S, where the states are labeled as \ng np ns)- The 
state 1 000) does not evolve under the Hamiltonian H of 
the considered three-level atoms (H |000) = 0), but only 
due to damping, and the Hamiltonian part of the evolu- 
tion takes place within the subspace spanned by 1 100) , 
1 010), and |001). Realistic atoms may require consider- 
ing more states, but in our problem the calculation of 
this 3x3 sub-matrix is equivalent to the calculation of 
£(t), P(t) and S(t) as carried out in this paper. Denoting 
{Ai} = {£, P, S}, the density matrix elements are simply 

n, = (AlA^t). (51) 

If the factorization condition (|4"9")l is satisfied, then 
r ij(to) = A*(to)Aj(to), which implies p(io) is a projector, 
and therefore that the initial density matrix describes a 
pure state. The other implication stating that if the ini- 
tial state is pure then Eq. (|4"9")l is satisfied follows directly 
from the definition (|51[) . Even more importantly, due to 
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Eq. ([50)) . the factorization is preserved at all times which 
implies that the density matrix takes the form 

p(t) = r(t)|000><000| + \m) m)\ , (52) 

where \ip(t)) = (£(t), P(t), S(t)) T is a (non-normalized) 
pure state. Therefore (up to a normalizing factor) the ini- 
tial state |001) evolves into the superposition £{t) 1 100) + 
P{t) |010) + S(i) |001). We see that the component of the 
single-photon state is exactly £{t). 



Appendix B 

We want to confirm Eq. (|18p which asserts that, for 
sufficiently large pulsewidths r, the quantum field £(t) 
simply follows P(t) 



£(t) = iwP(t-6t)/K 



(53) 



with a constant time delay St = 1/k. This relation is 
very similar to the relation between £{t) and P(t) in the 
bad cavity limit k w (which is not fulfilled in our 
case) except for the delay St. This similarity is related 
to the relation (|17j) in the frequency domain. In the time 
domain, we will show that it is closely related to the fact 
that the single photon escapes the cavity in a very early 
stage of the retrieval sequence. 

Let us assume that the read pulses are parameterized 
as Q(t) = Oo/(t/r). By using the dimensionless quanti- 
ties t/r, fcr, . . ., the number of parameters in Eqs. (fTj)- 
^ is effectively reduced by one, which makes the scaled 
quantities more natural. We can interpret Eq. (fl"8|) as 
the definition of a new function St(t). If St(t)/r is small 
at all times, then the conclusion St = 1/k directly follows 
that from Eqs. ((1) and (|T5|) . Therefore, we only need to 
show that 5t(t)/r is indeed small at all times for suffi- 
ciently long read pulses. This can be explicitly shown 



using the solution (|30)) . The implication that 1/k be- 
comes a short time scale in this limit also indicates that 
the single photon has been emitted very early in terms 
of relative times u = t/r. 

For long pulses, the real part of the argument of the 
exponential function in Eq. (f3T)| acquires very large nega- 
tive values which indicates that the photon emission has 
been terminated. The time interval during which the 
photon has been emitted can be estimated from the con- 
dition 



Re 



du' 



(54) 



Therefore, increasing r at fixed tto results in a photon 
emission at shorter relative times u = t/r. In the limit 
of very long pulses, all terms ~ fi/fi are negligible. Con- 
sequently, according to Eq. (|29|) . 



1 d£ 
£~dt 



o 



(55) 



p 2 + \n\ 2 

holds at all times during the photon emission. Combining 
Eqs. ((15)1. (JTJ, and ((SB), we get that 



(£(u + St(t)/r) - £(u)) /£{u) -> 



(56) 



is satisfied for all t during the photon emission. There are 
two logical possibilities consistent with the last equation: 
either St(t)/r is small at all times or St(i) relates distant 
instants with the same \£\ and phase 9g of the quantum 
field. The latter possibility is very restrictive and would 
be essentially an additional condition to the equations of 
motion ([I])-© since such distant correlations cannot nat- 
urally come from strictly local propagation equations. In 
fact, all panels in Fig. 2)clearly show that the points with 
the same \£\ have different phases so that the only con- 
sistent possibility derived from Eq. (|56[) is that St{t)/r 
is small for all t and consequently St = 1/k. 
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